NASA TECHNICAL NOTE 


i 




MSA 

TIT 

D-79UO 

C © 1 



NASA TN T 


<CdpYr re 
KlRTLA^D_ AJr g 


fc-* 

LU 

ID 

JZ 

DJ 

Ln 





UO 


'i- 

~A /FORMAL STRUCTURE FOR ADVANCED 

l f 

AUTOMATIC FLIGHT-CONTROL SYSTEMS] 

George Meyer and Luigi Cicolani 

Ames Research Center 
Moffett Field , Calif. 94035 

3* 

^NATIONAL AERONAUTICS AND SPACE ADMINJs¥R*#CTN • WASHINGTON, D. C. • 


I 




TECH LIBRARY KAFB, NM 



1. Report No. 2. Government Accession No. 

TN D-7940 

□ 1334*15 

4. Title and Subtitle 

A FORMAL STRUCTURE FOR ADVANCED AUTOMATIC FLIGHT-CONTROL SYSTEMS 

5. Report Date 

May 1975 

6. Performing Organization Code 

7. Author(s) 

George Meyer and Luigi CicOlani 

8. Performing Organization Report No. 
A-5710 

10. Work Unit No. 
501-03-11 

9. Performing Organization Name and Address 

Ames Research Center 
Moffett Field, Calif., 94035 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

An effort is underway at Ames Research Center to develop techniques for the unified design 
of multimode, variable authority automatic flight-control systems for powered-lift STOL and 
VTOL aircraft. This report describes a structure for such systems which has been developed to 
deal with the strong nonlinearities inherent in this class of aircraft, to admit automatic 
coupling with advanced air traffic control requiring accurate execution of complex trajectories, 
and to admit a variety of active control tasks. The specific case being considered is the 
augmentor wing jet STOL research aircraft. 


17. Key Words {Suggested by Author(s)) 

Handling qualities 
Flight controls 
Autopilot 


18. Distribution Statement 

Unclassified - unlimited 

STAR Category - 08 


19. Security Classif. {of this report) 

20. Security Classif. {of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

39 

$3.75 


For sale by the National Technical Information Service, Springfield, Virginia 22151 



TABLE OF CONTENTS 


Page 

SYMBOLS v 

SUMMARY 1 

INTRODUCTION 1 

BASIC COMMANDS TO AUTOMATIC FLIGHT-CONTROL SYSTEM 4 

TRACKING ACCURACY 7 

EQUATIONS OF MOTION 8 

AUGMENTOR WING JET STOL RESEARCH AIRCRAFT 13 

The Trimmap 16 

Perturbation Controller 18 

Angular Acceleration Controller 23 

Trajectory Command Generator 27 

PROPOSED STRUCTURE FOR ADVANCED AUTOMATIC FLIGHT-CONTROL SYSTEMS 31 

CONCLUSIONS 32 

REFERENCES 33 


iii 




SYMBOLS 


A direction cosine matrix, actual attitude of the aircraft with respect 

03 to inertial space 

A direction cosine matrix, commanded attitude of the aircraft with 

respect to inertial space 

A direction cosine matrix, commanded velocity axes with respect to 

inertial space 

b wing span 

a wing chord 

Cp drag coefficient 

commanded drag coefficient 

Cj cold thrust coefficient 

lift coefficient 

C T commanded lift coefficient 

Lo 

C^ a moment coefficient vector with respect to body axes 

^Mao commanded moment coefficient with respect to body axes 

C total force vector coefficient 

s 

D drag 

E * (<f>) elementary rotation about axis i through angle (f> 
f total aerodynamic and propulsive force in inertial coordinates 

f{x>u) right-hand side of system state equation 
g acceleration of gravity 

g(x,x) trimmap 

h a body coordinates of total angular momentum 

h{e) right-hand side of transition dynamics 
K(x y x ) feedback gain schedule 
L lift 
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m aircraft mass 

M body coordinates of aerodynamic and propulsive moment 

Q dynamic pressure 

R inertial coordinates of position vector 

R* inertial coordinates of position vector commanded by air traffic 

control 

R inertial coordinates of position given by command generator 

S w wing area 

t time variable 

ST throttle 

T cold thrust 

o 

T ^ hot thrust 

u control vector 

u a body coordinates of unit vector along relative velocity vector 

u inertial coordinates of unit vector along relative velocity vector 

s 

V airspeed (true airspeed) 

body coordinates of relative velocity vector 
v • measured body coordinates of relative velocity vector 

y g inertial coordinates of relative velocity vector 

7 g inertial coordinates of aircraft velocity vector 

7 * inertial coordinates of velocity commanded by air traffic control 

V qq inertial coordinates of velocity commanded by command generator 

V inertial coordinates of aircraft acceleration vector 

s 

V* inertial coordinates of acceleration commanded by air traffic control 

s 1 

7 ge inertial coordinates of acceleration commanded by command generator 

Ki inertial coordinates of acceleration input to trimmap 


vi 



V inertial coordinates of acceleration modifications due to perturbation 

controller 

Wg inertial coordinates of wind 

w s inertial coordinates of estimated wind 

x system state 

a angle of attack 

3 angle of sideslip 

glide-slope angle of relative velocity vector 

<5 elevator command 

eo 

& F flap angle 

<5. column matrix with 1 in ith row and 0 in the other two rows 

6 rudder command 

TO 

6^ throttle command 

<5 wheel command 

wo 

n variables of unsteady aerodynamics 

0 pitch angle 

v nozzle angle 

P density of air 

a side-force angle 

( f> roll angle 

ijj yaw angle 

body coordinates of aircraft angular velocity 
co^ bandwidth 
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A FORMAL STRUCTURE FOR ADVANCED AUTOMATIC FLIGHT-CONTROL SYSTEMS 


George Meyer and Luigi Cicolani 
Ames Research Center 


SUMMARY 


An effort is underway at Ames Research Center to develop techniques for 
the unified design of multimode, variable authority automatic flight -control 
systems for powered-lift STOL and VTOL aircraft. This report describes a 
structure for such systems which has been developed to deal with the strong 
nonlinearities inherent in this class of aircraft, to admit automatic coupling 
with advanced air traffic control requiring accurate execution of complex 
trajectories, and to admit a variety of active control tasks. The specific 
case being considered is the augmentor wing jet STOL research aircraft. 


INTRODUCTION 


Government and industry are investing substantial resources in developing 
new aircraft configurations required to meet the needs of the nation in the 
1980 f s and beyond. Present indications are that automatic flight-control sys- 
tems will play a significant role in this development. The basis for such a 
forecast is a combination of three factors. 

1. The mix of aircraft types such as VTOL, STOL, CTOL, and SST will 
require an advanced air traffic control (ATC) system. The accommodation of 
many aircraft covering a wide spectrum of speeds and maneuverability and at 
the same time satisfying stringent environmental constraints can be achieved 
only if the ATC has at its disposal a sufficiently large set of trajectories. 
Accurate execution of any one of a large set of complex trajectories will 
require a powerful automatic flight-control system that uses the maximum capa- 
bility of each aircraft type. 

2. Current work aimed at providing aircraft for short-haul transportation 
is developing the powered- lift technology. Among the concepts being considered 
are the augmentor wing, tilt rotor, lift fan, and externally blown flap. In 
all cases, the wide range of lift coefficient is achieved by inflight modifi- 
cations of the aircraft configuration. 

These modifications result in drastic changes in control characteristics 
of the aircraft; particularly in the high lift transition and landing con- 
figurations, the aircraft response to control inputs is very nonlinear. More- 
over, the presence of powered- and direct-lift generators increases the total 
number of controls available to the pilot who must continually make decisions 



on control techniques. Accurate, unaided manual tracking of complex trajec- 
tories by manipulating a large set of interacting controls of an aircraft 
whose control characteristics are nonlinear and rapidly changing represents an 
unacceptably high pilot workload. Automatic flight control can reduce the 
pilot workload to an acceptable level by integrating control functions so as 
to generate desirable handling qualities without reducing the performance of 
the aircraft as an element of the advanced civil air transportation system. 

The advantages of automatic flight control are potentially even more substan- 
tial in military applications of STOL and VTOL aircraft. Both the advanced 
military STOL and the Sea Control Fighter VTOL must utilize to the fullest the 
maneuvering capability of the basic aircraft. The tracking of complex trajec- 
tories must be sufficiently accurate to properly execute a mission, and the 
pilot workload associated with flying must not adversely affect his ability to 
perform other tasks. Again, the maneuverability, accuracy, flexibility, and 
level of pilot workload can be improved with automatic flight control. 

3. The rapidly advancing technology of sensors, actuators, and electronic 
components is approaching the point when servomechanisms with reliability com- 
parable to that of a wing can be built and maintained economically. Conse- 
quently, the conventional direct mechanical systems composed of cables, push 
rods, bell cranks, and mixers that link the pilot to control surfaces can be 
replaced by fly-by-wire systems. Although fly-by-wire technology itself 
offers several advantages over the conventional mechanical control systems, 
the real goal lies in the application of active control technology (ACT) to 
future aircraft. The key idea of ACT is the integration of control with aero- 
dynamics, structures, and propulsion early in the design cycle of the aircraft. 
Studies have shown that significant reductions in induced drag and structural 
weight, improvements in passenger comfort, and reduction of flight hazards can 
be achieved with ACT. These benefits are possible due to (a) a reduction in 
the sizes of stabilizing surfaces, with stability provided by dynamically 
controlling movable surfaces rather than statically with larged fixed surfaces 
as in the conventional designs; (b) reductions in structural strength require- 
ments by applying maneuver load alleviation and gust load alleviation; 

(c) improvement of ride qualities by a ride quality control system; and 

(d) reduction in the occurrence of inadvertent flight hazard through automatic 
limitation of flight conditions. These and other ACT concepts are currently 
being developed. A total automatic flight-control system is required to inte- 
grate all these control functions with the autopilot. 

Thus, indications are that automatic flight-control systems will play a 
significant role in the development of future aircraft. Of course, these 
systems were needed in the past, but the designer was severely limited by the 
characteristics of available transducers and, particularly, by the small 
inflight computational capacity. However, rapid advances have resulted in a 
large variety of accurate and reliable devices, while the capacity of digital 
flight computers has increased phenomenally and continues to increase. As a 
result, the designer is now limited primarily by the available methodology 
for the design of automatic flight-control systems. 

The most severe limitation of the existing design techniques is their 
extreme reliance on linear perturbation models of the aircraft. So long as 
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nonlinear effects are of minor significance, these techniques are quite ade- 
quate. But as nonlinearities become prominent because of either increased 
system accuracy requirements or the physics of force and moment generation in 
the powered-lift configurations, linear methods become less tractable. Many 
perturbation models are needed to cover the flight envelope adequately. Even 
the procedure for choosing reference trajectories about which to perturb is 
unclear at present, and controls corresponding to these trajectories that trim 
the aircraft cannot be generated easily or accurately by means of perturbation 
techniques. Logic must be provided in the flight computer for switching the 
perturbation control gains and reference controls as the aircraft leaves the 
domain of validity of one perturbation model and enters another. The result 
is a design that is complex in concept and implementation so that analyses of 
closed-loop sensitivity to modeling errors and subsystem failures are exceed- 
ingly difficult and not very convincing. 

Design techniques are needed of sufficient generality to be applicable to 
a large set of aircraft types with nonlinear dynamics and multiple redundant 
controls. The techniques must admit an effective tradeoff between tracking 
accuracy requirements on the one hand and requirements imposed on the capacity 
of the flight computer and on the a priori knowledge of system dynamics on the 
other hand. The techniques must be nearly algorithmic to permit tradeoff 
studies early in the aircraft design cycle when many alternative aircraft con- 
figurations are being considered. Techniques are needed for integrating a 
variety of active control functions with an autopilot having a multitude of 
modes and for coupling the autopilot automatically with the air traffic con- 
trol. Finally, these design techniques must result in designs sufficiently 
simple to admit an effective reliability analysis. 

An effort is underway at Ames to develop the methodology for the design 
of advanced flight-control systems. This report describes the progress made 
in the first segment of this program, namely, the formulation of an overall 
logical structure for multimode, variable authority, automatic flight-control 
systems. The proposed structure consists of five major subsystems: (1) The 

force trimmap trims the aircraft to any admissible time-varying acceleration 
vector. One of the outputs of the force trimmap is the possibly time-varying 
trim attitude. (2] The attitude control system generates commands to the 
moment-generating control surfaces and thereby forces the aircraft attitude to 
follow the input from the force trimmap. (3) The wind estimator provides 
estimates of the aircraft velocity vector relative to the air mass which are 
needed in the force trimmap and attitude control system calculations. (4) The 
trajectory perturbation controller closes the loop around the inaccuracies of 
the force trimmap, attitude control system, and wind estimator. The result is 
a trajectory acceleration vector controller whose input-output relation between 
the commanded acceleration and actual aircraft acceleration is essentially an 
identity, provided the input is flyable and its bandwidth is suitably 
restricted. (5) The trajectory command generator transforms the inputs from 
the air traffic control or the pilot into trajectories whose acceleration is 
consistent with the limitations of the trajectory acceleration controller. 

The basis for the proposed structure as well as its feasibility, benefits, and 
limitations are discussed. The internal structure of the five major subsys- 
tems is presented in some detail to clarify the intent of each subsystem. 
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The augmentor wing jet STOL research aircraft is used as a specific example. 

It is emphasized, however, that the objective of this report is not to present 
a complete automatic flight-control system for a particular aircraft, but 
rather to propose an overall logical structure for such systems. 

BASIC COMMANDS TO AUTOMATIC FLIGHT-CONTROL SYSTEM 


The boundary of the system considered here is shown schematically in 
figure 1. In the following discussion, the automatic flight-control system is 

the complete control system of the air- 
craft. It consists of all sensors, 
actuators, and control logic. The set 
of sensors measures aircraft motion and 
includes devices that are onboard as 
well as ground-based systems such as 
Figure 1.- Elements of automatic the MLS (when available). The function 
flight control system. of the control logic is to operate on 

the data from the sensors and commands 
from the (generalized) air traffic control (ATC) and thereby to generate com- 
mands to the actuators which, in turn, control the aircraft. The degree of 
automation of the control logic ranges from the fully automatic mode, in which 
the actuators are completely under the control of the flight computer, to the 
fully manual mode, in which the actuators are controlled exclusively by the 
pilot. Between these extremes, there is a spectrum of modes with specific 
functions such as handling quality control, ride quality control, gust load 
alleviation, maneuver load control, and a variety of autopilot modes such as 
autothrottle, altitude hold, heading capture, etc. Of course, combinations of 
such elementary modes may also be required. In addition, the control logic 
must be able to detect failures in various subsystems and switch (when neces- 
sary) to the next safest control strategy. The subject of this report is the 
design of such control logic. (The estimation problem associated with sensors 
and the design of fly-by-wire systems is not discussed.) 

The basic input to the control logic is the trajectory to be followed by 
the aircraft. The trajectory may be commanded explicitly by the ATC or 
implicitly by the pilot. A simple case, conceptually, occurs when (based on 
wind estimates, the capabilities of the aircraft, and other considerations) 

ATC selects a flyable trajectory to be followed by the aircraft. Generally, 
the set of admissible trajectories consists of a sequence of continuous seg- 
ments defined on relatively long (e.g., greater than 10 sec) intervals of time 
(ref. 1). Often the segments belong to a small set (e.g., lines and circles), 
in which case only the parameters and duration of the segments are transmitted 
to the aircraft and the commanded trajectory is reconstructed onboard. Other- 
wise, the coordinates of the trajectory are transmitted to the aircraft con- 
tinuously. In either case, the segments are defined on intervals of time; 
hence, position, velocity, and acceleration vectors corresponding to the com- 
manded trajectory are available to the control logic continuously. Moreover, 
since the motion of the aircraft in inertial space (a flat nonrotating earth 
is assumed throughout for simplicity) is of prime concern, inertial 



4 






coordinates of these vectors are considered as fundamental. The situation is 
essentially the same whenever the aircraft is commanded to coincide with a 
moving target as, for example, a carrier landing or docking with another air- 
craft or a missile intercepting another object. 

The pilot is an alternative source of commands. Of course, if he feeds 
the trajectory parameters into the autopilot either as a voice command from 
ATC or on his own initiative, he may be considered part of the ATC. However, 
many of the commonly used autopilot modes such as heading hold, altitude hold, 
autothrottle, glide-slope capture, control wheel steering, etc., generate the 
commanded trajectory only implicitly and often incompletely. Nevertheless, in 
most cases, an appropriate equivalent ATC trajectory can be constructed to 
represent the pilot command. The trajectory may contain a number of free 
parameters which the control logic can be instructed to ignore. Consequently, 
most commands concerning the motion of the aircraft center of mass may be con- 
sidered, at least conceptually if not in actual mechanization, to be generated 
in a standard form by the generalized ATC. 

In view of the preceding discussion, the following decision is made con- 
cerning the structure of the automatic flight-control system: 

Decision 1: The basic command to the automatic flight-control system is 

a concatenation of continuous segments T^, each of which is given by 



where the 9-tuple consists of 
inertial coordinates of commanded 
inertial position (7?^) , velocity 
(V s ) , and acceleration (7^) 
vectors . 

The complete trajectory may 
have discontinuities across the 
boundaries of the intervals Tfc. For 
example, all coordinates are discon- 
tinuous at t = in figure 2. The 
control logic must synthesize a 
transition trajectory consistent 
with the limitations imposed by 
dynamics. Another possibility is 



r* = ... r£ r £ +1 ... with 
total discontinuity at t - 
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Figure 3.- ATC command, T 
to land. 



illustrated in figure 3, which represents 
the vertical channel of the command to 
land. The segment r£ correspond to 
altitude variation while the aircraft is 
on the glide slope. The flare initiation 
altitude occurs at t = 0, at which time 
the segment r £ +1 is commanded. Thus 
there is a discontinuity in both position 
and velocity at t = 0. Again, the con- 
trol logic must synthesize an appropriate 
transition (flare) trajectory. 

As already noted, some parameters of 
P* the commanded trajectory may be free. In 
^ +1 * particular, all nine coordinates need not 
be always tracked. For example, consider 
the three segments shown in figure 4. 
Segment rj£ represents the command to 
track a four-dimensional trajectory with 


constant altitude h . 


At t = t, , the 


Figure 4.- Encounter with heavy 
turbulence. 


aircraft encounters Seavy turbulence. 
Depending on the severity of the turbu- 
lence relative to the limits imposed on 
aircraft dynamics, tracking may have to 
be relaxed from position, to velocity, or 
acceleration and, in the extreme case, 
only the attitude of the aircraft will be 
tracked while the ATC command is ignored 
completely. As a result, the aircraft is 


allowed to drift along some trajectory r£ +i away from the true ATC command. 
As turbulence subsides at t = t 2 , four -dimensional tracking can be resumed. 
However, because of the errors accumulated in the interval there will 

be, in general, a discontinuity in all coordinates of the command at t = t 2 . 
The control logic must synthesize an appropriate transition trajectory to 
bring the aircraft back on r £ +2 * The situation is essentially the same when 
the set of operating sensors changes with time or when the constraints imposed 
on the aircraft dynamics change perhaps because of failures in certain sub- 
systems. 

Based on the preceding discussion, the following decision is made con- 
cerning the formal structure of the control system. 

Decision 2: The control logic contains a command generator that synthe- 

sizes trajectories 


r = l 


= 


r o (t)\ 

SO K ' 


V ( t ) 

SO K J 


T sc ^i 


vt 
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which are flyable at all times by 
the aircraft with the available 
sensors and actuators and with 
existing constraints imposed on 
dynamics . 

For example, the output of 
the command generator corresponding 
to the case in figure 3 is shown in 
figure 5. Note that there is no 
discontinuity in r at t = 0; while 
there is a discontinuity in r* at 
the same instant. (The command 
generator is discussed further 
later in the report.) 



Figure 5.- Output of the command generator 
for a landing maneuver. 


TRACKING ACCURACY 


Clearly, one of the essential attributes of a control system is that it 
be as simple as possible, both in concept and in mechanization. The level of 
complexity, however, is determined ultimately by accuracy requirements. If 
the required accuracy is low, then many details of the aircraft equation of 
motion may be suppressed, and a simple model usually leads to a simple control 
law. As accuracy requirements are increased, a more detailed representation 
of aircraft dynamics becomes necessary. The model increases in complexity. 

The dimension of its state space increases as more dynamical elements are 
accounted for. New cross -coup ling links appear. The number of parameters 
increases with finer representation of nonlinearities. All this increases the 
complexity of the control system. Moreover, the design methodology may have 
to be changed completely with an increase in accuracy requirement. As a 
result, tradeoff studies may become intractable. However, such studies are 
essential in the design of automatic flight-control systems advanced aircraft. 
Of particular interest is the tradeoff between the required capacity of the 
onboard computer and trajectory tracking accuracy. Hence, a single design 
methodology must be developed in which tracking accuracy is a variable. 

The natural evolution of an AFCS for a new, possibly experimental, air- 
craft is by means of a sequence of refinements. For safety, initial flight 
tests are restricted to relatively simple maneuvers and to correspondingly 
simple modes of the control system with minimal authority and tracking accu- 
racy. As flight data accumulate and good estimates of critical aircraft 
parameters become available, the set of maneuvers is expanded until, finally, 
it coincides with the designed flight envelope. Thus, the automatic flight- 
control system must be designed to allow a spectrum of tracking and modeling 
accuracies. 

A spectrum of accuracy, rather than a single level, is also needed for 
normal aircraft operation. For example, in cruise, altitude tracking is 
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obviously not as significant as during a landing maneuver and can be traded 
for, say, ride quality. 

Therefore, the following decision is made concerning the structure of the 
control system. 

Decision 3: The tracking accuracy enters the control logic as an inde- 

pendent variable, both during design as well as in normal operation. 

The accuracy of a control system is ultimately limited by the accuracy of 
the navigation system. Hence the accuracy of the latter serves as an estimate 
of an upper bound on the former. The RAINPAL system (ref. 2) is one of the 
most accurate, flight- tested, navigation systems currently available. A com- 
parison of RAINPAL errors with allowable errors for CTOL and SSV is given in 
table 1. In the remainder of this report, the RAINPAL errors are taken as the 
upper limit on the trajectory tracking requirements. 


TABLE 1.- COMPARISON OF RAINPAL NAVIGATION ERRORS WITH 
ALLOWABLE ERRORS FOR CTOL AND SSV 


Component 

RAINPAL navigation 
error standard 
deviations 

Navigation errors 
allowable for CTOL 
automatic landing 
systems 

Navigation errors 
allowable for the 
SSV autoland 
system 

X 

0.9 ±0.6 m (3 ±2 ft) 

132 m (433 ft) 

43.2 m (139 ft) 

Y 

1.2 ±0.6 m (4 ±2 ft) 

2.38 m (7.79 ft) 

1.51 m (4.97 ft) 

Z 

0.9 ±0.6 m (3 ±2 ft) 


2.46 m (8.08 ft) 

• 

X 

0.15 ±0.06 m/sec 
(0.5 ±0.2 ft/sec) 


1.76 m/sec 

(5.77 ft/sec) 

• 

Y 

0.3 ±0.15 m/sec 
(1 ±0.5 ft/sec) 


0.88 m/sec 

(2.89 ft/sec) 

• 

Z 

0.15 ±0.06 m/sec 
(0.5 ±0.2 ft/sec) 


0.088 m/sec 

(0.289 ft/sec) 


EQUATIONS OF MOTION 


Let the inertial coordinates of the aircraft position and velocity with 
respect to the runway axes (flat, nonrotating earth is assumed throughout) be 
denoted by R s and V s , respectively. Then 
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R 

s 


V 

s 


(3) 


where (*) denotes differentiation with respect to time t . Aerodynamic forces 
and moments depend on the velocity of the aircraft relative to the air mass. 
Let W s denote the inertial coordinates of the wind velocity. Generally, W s 
is a complicated function of position and time which may vary significantly 
over the dimensions of the aircraft. Let 


W s = W 3 + 6w e (i? s + p s>V 


(4) 


where the first and second terms consist of wavelengths longer and shorter, 
respectively, than the aircraft dimensions, and r s is position-referenced to 
the aircraft center of mass. The inertial coordinates v Q of the aircraft 
velocity relative to the air mass are defined in this report by 


V 

s 


- ^ 


(5) 


where w s = w s (R s ,t) is the average wind at the aircraft center of mass. Wind 
shear across the aircraft is ignored here. Polar coordinates of relative 
velocity are defined in a standard manner according to figure 6. Thus 


y s = vu s (63 

and 

where (f>) is an Euler rotation about 
axis i y ( is the transpose of ( ) , 
and is the column with 1 in the ith 
place and 0 in the other two. In the 
absence of wind, is the aircraft 
heading angle and y v is the glide-slope 
angle. 

In che aircraft body axes, the 
relative velocity is given by 


where 



Figure 6.- Definition of airspeed, 
v , heading angle \p v and glide 
slope angle 


u a = E 2 C a)^(6)6 1 


(9) 
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Conversely, 
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a = tan -1 [w a (3)/M a (l)] 
3 = sin ' 1 [u (2) ] 


( 10 ) 


where a and g are the angle of attack and sideslip angle as normally defined 
(ref. 3) . 

The attitude of the aircraft body axes with respect to the runway axes 
are defined by the direction cosine matrix A as . If Euler angles are used in 
the 3-2-1 sequence, then 


A as - V»Ve)£ 3 C« 


( 11 ) 


The attitude can also be defined in terms of the angles associated with the 
direction of the relative velocity vector as 


as 


= E. 


(a)£( 3) 


(OS, CYJS, OU 


V 


( 12 ) 


where <f> y is the angle of roll about the relative velocity vector v s . A block 
diagram representation of equation (12) is given in figure 7 for future 

reference. The term "heading" 

Velocity Wind tunnel Stability Body 


Axes: Runway Heading 


I- 


E 3 (*,)|_*|E 2 (y v > 


E|(*„) 


refers to the heading of the rela- 
tive velocity vector. 


Subscript s 


h v wt st a 


Figure 7.- Block diagram representation 
of equation (13) . 


Let the body coordinates of 
aircraft angular velocity with 
respect to the runway be denoted by 
co a . Then (see, e.g., ref. 4), 



= SfrJA 


as 


where. 


for any column x = 


(#1 ,x 2 




T 


S(x) = 



s c. 


0 


-X, 



(13) 


(14) 


Equations (3) and (13) are the kinematic components of the aircraft equations 
of motion. 


Let the inertial coordinates of total aerodynamic and propulsive force be 
denoted by f s , and let m and g be the aircraft mass and acceleration of 
gravity, respectively. Then 
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r 


V 

s 



+ & S 3 


( 15 ) 


The total aerodynamic and propulsive force is most directly expressed in terms 
of coordinates with respect to the wind-tunnel axes. Thus the total force 
along the relative velocity vector, henceforth to be called total drag, 

°wt ■ -«W l < 16 > 

where the dynamic pressure 

Q = \ pv 2 (17) 


for which p is the density of air, is the wing area, and Cp is the total 
drag coefficient. The total force perpendicular to the relative velocity 
vector, henceforth to be called total lift. 





(18) 


where C p is the total lift coefficient and a is defined in figure 8. Note 
that the present definition of the total lift coefficient includes the side 
force, and both total lift and drag 
coefficients include the effects of 
thrust. Generally, 


C D = 

Cj^ = C L {a 9 $yU-yT\) 

o = a(a,3,u;n) 

where u represents the available con- 
trols such as flaps, throttle, ele- 
vator, rudder, ailerons, etc., and n 
represents the dynamic variables such 
as a, 3, u> a , etc., and other variables 
such as air temperature and density. 

The inertial coordinates of the 
total aerodynamic and propulsive force 
are given by 



L 


■*■2 


3 



f Q ~ QS w C q (20) Figure 8.- Definition of lift 

vector used in report. 

where the total force vector coeffi- 
cient is 
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( 21 ) 


C S = “ il L®2 c °° s 3 cwc ‘ 7 i> 6 1 + C lfl Ca)6 s] 

The dynamic equation for rotation is given by 

J) = J~ l [M + S(i» )h ] 
a a L a v a a i 

where J a is the aircraft moment of inertia in body axes, M a is the total aero- 
dynamic and propulsive moment, and h a is the total aircraft angular momentum. 
When the angular momentum of spinning parts is negligible. 


h 


a 


J u 
a a 


(23) 


The moment vector is defined in terms of the moment coefficients in the usual 
manner : 


M 


a 



0 

a 

0 



(24) 


where b and c are the span and mean chord, respectively, of the wing. 
Generally, 


C., = C 
Ma m 



Figure 9.- Main information flow in 
force and moment generation. 


,g,u;ri) (25) 


The data represented by equa- 
tions (19) and (25) are considered as 
the fundamental source of information 
concerning the total aerodynamic and 
propulsive force and moment. In the 
remainder of this report, these data 
are assumed to be given to various 
levels of accuracy. The information 
flow involved in force and moment 
generation is shown in figure 9. 


Equations (3), (13), (15), and (22) are the fundamental components of the 
system state equation. Other effects such as the dynamics of actuators and 
sensors may be adjoined to these equations as the need arises to obtain the 
complete state equation modeling the aircraft. 
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AUGMENTOR WING JET STOL RESEARCH AIRCRAFT 


A specific aircraft is described here to motivate and aid in the follow- 
ing discussion. Note that, although the discussion in the remainder of the 
report is directed toward this specific aircraft, the essential concepts are 
applicable to other types of aircraft. 


The augmentor wing jet STOL 
research aircraft (AWJSRA) is a 
de Havilland C-8A "Buffalo” modified 
according to the general arrangement 
shown in figure 10. The wing area 
S w is 80.36 m 2 (865 ft 2 ) and the 
maximum gross weight is 20,400 kg 
(45,000 lb). The aircraft is 
powered by two turbofan engines. 

The relatively cold flow from the 
front fans is ducted through the wing 
and fuselage to the augmented jet 
flap and blown ailerons. The 
arrangement of the augmentor flap is 
shown in figure 11. The entire flap 
unit pivots about the main hinge 
point. No provision is made in this 
installation to retract the upper 
flap units into the main wing con- 
tour. The Coanda surface serves to 
deflect the (cold) flow from the 
nozzle. The augmentor chokes at the 
trailing edges of the main flaps con- 
trol the lift generated by the flaps. 
The two outboard flap chokes are 
used to control roll and all four 
chokes are used to spoil lift after 
landing. 

The hot gas from the two turbo- 
fan engines flows through two pairs 
of nozzles that can be rotated in 
flight to provide vectoring of the 
hot thrust through a range of 98°. 

All nozzles are connected to move in 
unison in response to a single nozzle 
angle command. The geometry asso- 
ciated with the hot thrust is shown 



Figure 10.- General arrange- 
ment of the modified C-8A. 



Figure 11.- Arrangement of the 
augmentor flap. 
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in figure 12. Since the aircraft center of 
gravity is not on the axis of rotation of the 
nozzles, the hot thrust generates a moment 
that depends on the nozzle angle v. The 
servos that control the nozzles are quite 
fast, being limited to 90°/sec. The hot and 
cold thrusts depend on the engine speed. The 
speed of both engines is controlled in unison 
by a single throttle command, 6^. The asso- 
ciated servo system is relatively slow with a 
bandwidth of approximately 1 rad/sec. The 
cold flow has a pronounced effect on the wing- 
body polars of the aircraft as shown in 
figure 13, where the independent variables are 
flap angle 8p, angle of attack a, and cold 
thrust coefficient Cj = T 0 /QS W . 


Of particular significance for the 

Figure 12.- Geometry o ot design of automatic trajectory tracking sys- 

thrust in body coordinates. terns is the large variation in the basic aero- 

dynamic characteristics of the aircraft (evident in fig. 13). Certainly, there 
is quite a significant change between cruise configuration (flap =4.5 ) and 
landing configuration (flap = 65°). But present indications are that the non- 
linearity is significant even over much smaller regions. For example. 



Figure 13.- Wing-body polars (C^C^) (a,^^) . 
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figure 14 shows the total lift and drag 
coefficients, including the effects of the 
hot thrust for a constant flap, throttle, 
and airspeed which correspond to a typical 
landing configuration with angle of attack 
and nozzle angle as the independent 
variables . 

Point A 1 represents equilibrium flight 
at 65 knots along the -7.5° glide slope. 
Point A 2 represents level flight at the 
same speed. The derivatives of the total 
force vector coefficients are also shown 
at these two points. As the aircraft is 
maneuvered between A^ and A 2 , the change in 
these derivatives may adversely affect 
closed-loop stability. Of additional con- 
cern is the magnitude of steady-state error 
in altitude which will result if linearity 
is assumed in this region. Thus, for a 
maneuver that takes the aircraft from A x 
to A 0 , the feedforward based on the linear 


Flap = 65° 
Throttle = 22° 



Figure 14.- Total force 
coefficient. 


model at A x will trim the aircraft to 
point A 12 giving an error in the required 
lift coefficient, A C^/C^ = 4.5 percent. This error will be absorbed by the 
altitude feedback, resulting in an error h e = 0.045^/co^, where g is the accel- 
eration of gravity and is the bandwidth of the altitude loop. But, because 
of limitations imposed by the unsteady aerodynamics on this bandwidth, 
a ] n - rad/sec at 65 knots; hence h e ^ 1.8 m (5.8 ft). Similarly, a maneuver 

that takes the aircraft from level flight (point A 2 ) onto the -7.5° glide slope 
(point Aj) by means of the feedforward based on the linear model at A 2 will be 
trimmed to point A 21 , resulting in a magnitude of altitude error h e > 4.9 m 
(16 ft). Of course, this hang-off error can be removed by means of an inte- 
grator in the altitude error channel. Because of bandwidth limitation, such 
trim corrections will be too slow for many maneuvers. Consequently, when such 
errors cannot be tolerated, even the relatively small transition as between 
A 1 and A 2 must be considered nonlinear for the purposes of generating trim 
controls. So, as usual, the practical concept of linearity is intimately 
related to accuracy requirements. 


The trim problem is further complicated by the presence of redundant con- 
trols. Thus, as shown in figure 14, the total lift and drag coefficients 
required for steady flight at 65 knots along the -7.5° glide slope (point A 1 ) 
can be achieved by the trim condition (a, v, <5^, 6^0 = (2.5,92,22,65). However, 
the same coefficients can be achieved by other combinations of controls. This 
redundancy is shown in figure 15, where point A 1 corresponds to the particular 
solution Aj in figure 14. The question (the redundant control problem) is 
which combination (a, v, 6/^,6^) to use in any given situation. 


The preceding discussion of the characteristics of the AWJSRA implies 
that careful attention must be given to the problem of trimming the aircraft. 
Present indications are that other aircraft types that use powered lift and 
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Figure 15.- Controls for one 
value of total force 
coefficient. 


undergo rapid changes in configuration 
are similarly nonlinear and have redun- 
dant controls. Hence the following 
decision is made concerning the struc- 
ture of the automatic flight-control 
system. 

Decision 4: The control logic con- 

tains a section in which the control 
redundancy is resolved and trim controls 
are generated continuously. This sec- 
tion of the control logic is referred to 
as the trimmap. 


The Trimmap 


To solve the trim problem, one must, in effect, reverse some of the 
information flow shown in figure 9. Thus, given the relative velocity v s and 
density p, and the commanded (trim) total force vector f so and moment vector 
M ac , the problem is to find the required trim controls u a . From equation (20), 
it follows that 

C = (QS ) _1 f (26) 

so ^ w J J so v J 


Since drag is defined as the force component along -u s , and lift is defined as 
the total force perpendicular to u g . 


'Do 


T 

~ u s C so 


Lo 


^so C so 


- w 2 


(27) 


where u s is defined by equation (6) . 

Two cases arise in the computation of the commanded side-force angle o c 
(see fig. 8), according to whether or not the commanded (trim) attitude of 
the aircraft A* s is completely defined outside the trimmap. 

< & 

If A oS is completely defined, then the commanded angle of attack a e and 
side-slip angle are defined because the relative velocity vector v s is not 
subject to control within trimmap; when in trim, v a = A os v s and a and 6 are 
defined by equation (10) . Consequently, the wind-tunnel coordinates (see 
fig. 7) of the total aerodynamic and propulsive force vector coefficient 
required for trim are given by 

c vto • C28) 
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while equation (21) implies that 


a a = arctan[<7 ytc (2), ~ C wto m 


(29) 


so the problem is reduced to that of partially inverting the basic data of 
equations (19) and (25) : 


^(a.e.Mj h) - c Dq = ° 

^(a.e.wjn) - = 0 

a(a,B,w;n) - = 0 

C^(“»3,M;n) - C Maa = 0 


(30) 


If /4^ s is incompletely defined in that the trimmap is free to select the 
commanded angle of attack, then this angle and the controls must be chosen to 
satisfy equation (27) with constraints (30) , which define the side-force angle 
a a . Then the commanded attitude can be defined according to equation (12), 
namely. 


where 


(31) 

<f> y = ~o Q + arctan[C ye (2) , -C VQ { 3)] (32) 


and (see fig. 7) 


vc 


= Z 2 (Y v )E 3 (t,)C, 


V ' sc 


(33) 


The commanded attitude can also be defined without the explicit use of 
the angles (<j> y , as follows. The unit vector along the lift vector is 
given by 


■ V,a * C Dc U B>' C Le (34) 

where C-q c and Cj^ Q are defined by equation (27) . Let the matrix "us be the 
rotation defined by 

U US = K’ - S( - u s )uP s’ (35:) 

where S is the vector cross-product operator defined by equation (12) . The 
rotation U us takes the axes that are initially aligned with the runway (iner- 
tial) axes into the attitude in which the relative velocity is along 6^^ and 
the lift vector is along -6g. Hence the trim attitude is also given by 
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(36) 


4 ■ 

The main flow of information in the automatic trim logic when the commanded 
attitude is incompletely defined is shown in figure 16, The primary inputs are 

the input translational accelera- 
tion V s j and the input angular 
acceleration a )qj. The output is the 
control, u, and the required trim 
attitude, A QQ . The solution of 
equations (30) (called the trim- 
map) is the core of the automatic 
trim logic. Within the trimmap, 
control redundancy is resolved and 
control strategy is modified in 
case of component failures. The 
trimmap provides a natural setting 
for monitoring the proximity of the 
aircraft to its performance limits 
and for protecting the aircraft 
from exceeding its design limits, 
that is envelope limiting. Further- 
more, the primary purpose of the 
automatic trim logic is to provide 
a priori open- loop information to 
the overall automatic flight-control system and thereby relieve the perturba- 
tion controller whose feedback is intended to control the uncertainties of the 
process as well as relatively insignificant details that are known but ignored 
in the construction of the trimmap. Thus, it is also within trimmap that the 
major tradeoff between complexity and computer capacity on the one hand and 
accuracy of performance on the other takes place. 

The perturbation controller is discussed in the next section. However, 
note that the relative velocity vector v s is used in the trim logic. Since 
wind contributes significantly to the relative velocity, estimates of the wind 
must generally be computed. For this reason, the following decision regarding 
the structure of the automatic flight-control system is made. 

Decision 5: The control logic includes a wind filter that estimates the 

inertial coordinates, w Q , of the wind vector. (The wind filter is discussed 
further later in the report.) 


m 



Figure 16.- Main information flow in 
automatic trim logic. 


Perturbation Controller 

The perturbation controller provides closed-loop, feedback control over 
details of the physical process not accounted for in the open-loop, feed- 
forward, trim control either because they are not known a priori or because 
they have been purposely ignored to simplify the open-loop control. For the 
purposes of discussion, let the system state equation be 
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x - f(x,u) 


(37) 


where x and u are the state and control, respectively, of appropriate dimension 
and, in addition, the control is restricted to a set U that may depend on the 
state. A trajectory [x Q {t ) , t e T] is flyable if, for all t in T, there is a 
control u Q {t ) such that 

i 0 t *) - (38) 

The trim problem (as discussed previously) is to find a control u Q that satis- 
fies equation (38), given that the trim (nominal) trajectory x Q is flyable. 

The solution will be an inverse of the state equation (37), namely, a function 
( g yF') y which we call the trimmap, so that for all (a?,x) in F, 

g (cc , x) ] = x (39) 

The corresponding trim control is given by 

u o = 8&o* x o) ( 40 ) 

Usually, trim refers to cases with constant u Q . Here the concept is general- 
ized to include open-loop controls that vary with time. As noted previously, 
when controls are redundant, the state equation (37) alone does not suffice to 
define the trimmap (g,F), and additional conditions must be introduced to 
resolve the redundancy. 

The trim problem may be difficult to solve, but, evidently, its solution 
to the required accuracy is the essential first step in any design of auto- 
matic flight-control systems. The next step usually is to design feedback 
control systems based on perturbation models. Thus, given a flyable trajec- 
tory {x Q ,x 0 ) trimmed by u Q according to equation (40), the linear model (41) 
is obtained for the perturbations &x = x - x Q and 6u = u - u Q : 



where the partial derivatives are evaluated along the nominal trajectory. 
Then the application of the methods of linear control theory (ref. 5) yields 
the perturbation control law 


Su = K q 6x (42) 

Since the coefficients in equation (41) depend on the nominal trajectory, the 
process must be repeated for a sufficiently large number N of nominal trajec- 
tories (x Q ,x 0 ) in F until the flight # envelope F is covered adequately. The 
result is a scheduled gain matrix K{x 0> x 0 ') and the complete control law is 
given by 


(43) 


U = g(x Q ,x o ) + K(x o ,x o )(x - x o ) 
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The structure of the control 
system selected in the conventional 
perturbation controller design is 
shown in figure 17. The distur- 
bance is represented by e. 

The many advantages of the 
methods of quadratic optimization 
of linear systems are well publi- 
Figure 17.- Selected structure in cized (ref. 5). Indeed, these 

conventional design. methods are very powerful. But when 

the state equation (37) is highly 
nonlinear, the procedure for choosing the proper set { (x 0 i>Xoi) }N of nominal 
trajectories to adequately cover the flight envelope F with perturbation models 
so as to apply these methods to each of them is, at present, rather unclear. 

Then there is the problem of how to obtain envelope limiting. The perturba- 

tion control 6u (see fig. 17) may, in response to a disturbance, force the air- 
craft outside its design limits even though the nominal control u Q is main- 
tained by the trimmap well within these limits. If envelope limiting is 
achieved by limits SU on 6m, then questions concerning the stability of the 
resulting nonlinear perturbation controller must be resolved for each of the N 
perturbation models. Since these limits on u are likely to depend on (x 0 ,x 0 ') , 
6U(x 0 ,x 0 ) must be stored in the control logic in addition to the gain matrix 
K(x 0 ,x 0 ) and the trimmap g[x 0 ,x 0 ') . If N discrete cases of these variables are 
stored, the dynamics of switching from one case to the next as the aircraft 
state moves through their domains of validity must be designed. For example, 
hysteresis may have to be introduced to prevent jitter when the aircraft is 
required to fly near the boundary separating these regions. Switching 
dynamics further increase the storage and computation requirements of the 
flight computer. Lastly, there is a plethora of problems associated with 
reliability. Typically, does the system remain stable when a column or row of 
the gain matrix K is lost? — which corresponds to the loss of a sensor or a 
control, respectively. Such questions concerning the structure of the feed- 
back are, as yet, difficult within the framework of quadratic optimization. 

In this respect, the older classical design techniques based on sequential 
loop closures that result in a nesting (hierarchy) of subsystems with decreas- 
ing bandwidth are more effective for designing fail-safe systems. 

Because of these considerations, the following decisions are made con- 
cerning the formal structure of the automatic flight-control system. 

Decision 6: The feedback is closed through the automatic trim logic. 

Decision 7: The structure of the control logic is hierarchical. 

The information flow implied by decision 6 is shown in # figure 18. The 
feedback is through perturbation Sx Q in the trim condition x Q . Suppose that, 
initially, x = x Q . In the absence of any modeling errors, the control 

u = g(x Q ,x) (44) 
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will maintain x = x Q . The tracking 
will be perfect even if, at some point 
in time, x Q is perturbed to x Q + 6x a , 
provided that ( x Q + Sa;^,#) is in F. 

The corresponding control is 

u = g(x Q + &x 0 ,x) (45) 



Now suppose that initially x - x Q =£ 0, Figure 18.- Structure of proposed 
but that the error can be removed by perturbation controller, 

means of a flyable trajectory. Then 

there is an x Q + Sx Q that will take x into x Q by means of the control law 
given by equation (45). That is, the feedback for controlling the process 
uncertainties can be closed through the automatic trim logic as in figure 18 
rather than after the trim logic as in figure 17. One immediate consequence 
is that envelope limiting is done within the trimmap. The other consequence 
is that, for any admissible Sx Q , the perturbation model is given anywhere 
inside F by 


Sx = Sx + e (46) 

o 

where the magnitude of the error e depends on the accuracy of the automatic 
trim logic. Thus the emphasis is shifted from the N perturbation models 
required to cover F to the construction of flyable perturbations in the com- 
manded trajectories. The latter task is considerably simplified by decision 7. 


Consider the block diagram in figure 19, which represents the automatic 
flight-control system as viewed from one level in the hierarchy, namely, that 
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of the acceleration controller. The function of the acceleration controller 
is to accept commands from the command generator, which is one level higher, 
and transform them into commands to the flap, throttle, and nozzle servos, as 
well as to the attitude control system, all of which are one level lower. At 
the level of the acceleration controller, the servo systems are represented by 
relatively simple, possibly linear and low-order, input-output relations, which 
are treated as specifications to be met in the design of these subsystems. Of 
course, the subsystems may be quite complicated internally. For example, the 
attitude control system may have its own automatic trim logic and perturbation 
controller and may rely on simple input-output relations that describe the 
control surface servos, which are another step lower in the hierarchy. 


The major blocks of the acceleration controller are the trimmap, wind 
filter, and compensator. The inertial coordinates V g j of the input accelera- 
tion vector are transformed by the trimmap (fig. 16) into commanded flap, 
throttle, nozzle, and attitude. 


The wind filter computes smoothed inertial coordinates w s of aircraft 
velocity relative to the air mass from body-mounted air velocity v sensors 
and from the inertial velocity V s and attitude A as of the aircraft. Note that 
only the inertial coordinates of wind are filtered, while V s is unaffected. 
Hence, in the absence of wind and, of course, sensor errors, v s = 7 g . 

* • • 

The input-output relation, V s j -*■ V g , where the terms are the inertial 
coordinates of aircraft acceleration vector, is given by (see eq. (46)), 


V 

s 


V sl + e 


(47) 


where e depends on the inaccuracies of the trimmap and the wind filter, the 
presence of unsteady aerodynamics such as a effects , and on attitude and other 
subsystem dynamics. The purpose of the perturbation controller is to close 
the loop around such effects and thereby reduce e to a tolerable level. Iner- 
tial coordinates of position, velocity, and acceleration errors are trans- 
formed into approximately longitudinal, lateral, and normal errors by means of 
the direction cosine matrix A vs computed from the commanded inertial velocity 
V sc ', the errors are weighted by constant gain matrices , K 2 , and K 3 commen- 
surate with the acceleration capacities of the aircraft in these directions. 
The result is filtered to ensure compatibility with the attitude control sys- 
tem and other subsystem dynamics. The corrective acceleration is transformed 
back into inertial space and added to command V sc to give input F s j. In this 
way, feedback is closed around the process uncertainties, e, so that 


V 


V 


sc 


(48) 


is sufficiently accurate if the acceleration V sc commanded by the command 
generator is admissible, namely, if (F se ,F sc ) is flyable and the bandwidth of 
V sc is suitably restricted. Coriolis terms due to the time rate of change of 
A vc may be included in the perturbation controller if necessary using the 
techniques of the next section. 
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Angular Acceleration Controller 


The discussion thus far has been concerned with controlling the motion of 
aircraft's center of mass. The concepts that led to the structure of the 
(translational) acceleration controller shown in figure 19 are also applied to 
formulate the structure of the angular acceleration controller. The function 
of the angular acceleration controller is to accept commands from the attitude 
command generator and transform them into commands to the wheel, elevator, and 
rudder servos. The structure is again hierarchical. The attitude command 
generator, one step above the angular acceleration controller, accepts attitude 
requests from the translational control system and, based on simple input- 
output representation, generates rotational trajectories [A as (£) ,a) e (t) 
as input to the angular acceleration controller. The control surface servos, 
one step below the angular acceleration controller, are represented by rela- 
tively simple, input-output relations. 

The structure of the angular acceleration controller is shown in figure 20. 
The major blocks are the moment trimmap, wind estimator, and attitude pertur- 
bation controller. The body coordinates u) a j of the input angular acceleration 
vector are transformed by the moment trimmap into commanded wheel 6w e , ele- 
vator Se a , and rudder 



Figure 20.- Structure of the angular acceleration controller for the AWJSRA. 

The wind estimator provides estimates of the body coordinates v a of air- 
craft velocity relative to the air mass, which are needed in the moment trim- 
map calculations. The structure is very similar to that of the wind estimator 
in the translational control system. There may be differences in detail 
because of different bandwidth requirements. 















The input-output relation, to a ^ -*■ u> a , is given by 


<0 = 

a 


CO 


ai 


+ e 


( 49 ) 


where e depends on the inaccuracies of the moment trimmap and the wind esti- 
mator, the presence of unsteady aerodynamics such as a effects and wind shear 
across the aircraft, and on the dynamics of control surface servos. The pur- 
pose of the attitude perturbation controller is to close the loop around such 
effects and thereby reduce the error e to a tolerable level. Following 
reference 4, attitude error is defined by the direction cosine matrix 


T 

A = A A 
ao as os 


(50) 


that represents the aircraft attitude relative to the commanded attitude, 
time derivative is 


The 


But (see eq. (14)), 


• « T 7 • V 7 

A ~ A A + A A 
ao as os as os 


A = S(co )A 
as K a as 


and 


A 

os 


= S(oo e )4 


OS 


where co a represents (in body axes) the aircraft angular velocity relative to 
the runway, and co 0 represents (in the commanded body axes) the commanded angu- 
lar velocity relative to the runway. Hence 

A = 5(w ~)A -A 5(w ) 
ao a ao ao o J 

= SO );4 - A SO )A T A 

a ao ao o J ao ao 

= SO - S(A CO )A 

a J ao *■ ao o J ao 

= S(co - A (o )A 

a ao o J ao 


or, equivalently. 


A = SO )A 
ac ea ac 


(51) 


where the body coordinates of angular velocity error are given by 


co = co - A co 
ea a ao o 


(52) 
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Therefore, the time derivative is 


co =o) - iSfo) )4 a) - A o) 

ea a K ea J ao o ao c 


~ 0) + S fm -4 0) 

a ^oa JUJ ea ao o 


where u ea = 4 ae o) (3 represents (in the body axes) the commanded angular velocity. 
The identity S(x)y = - S (y^x was used in the last step above. Thus the body 
coordinates of angular acceleration are expressed in terms of the command and 
perturbations as 


03 =4 03 - S f 03 ) 03 +03 

a ao o v oa J ea ea 


(S3) 


Note that, since no small signal approximations are used to derive equa- 
tion (53), it is universally valid. 

Equation (53) can be interpreted as defining the required angular accel- 
eration of the aircraft (which is the input to the moment trimmap) , so 
that the command is executed with perturbation {A ao > ^ea > ^ea) • equation 

connecting w ea to (4 ac ,aj ea ) closes the loop around the perturbations. Thus 
the remaining problem is to synthesize a control law h = so that 

the system 


A = S( a) )A 
ac ea J ao 

(54) 


has an acceptable relaxation transient response, [A ac (0) ,w ea (0)] (T,0). 

This problem is treated in some detail in references 4 and 6. For example, 
the attitude error is defined by 



where (<££j) = A ac . According to Euler's theorem on rotations, every attitude 
can be attained with a single rotation. Let <f> be the angle of A ao , and let o 
be the unit vector along the axis of A ao . It can be shown that e = (sin <j>)c. 
Thus, for small attitude errors (<f> < 0.5 radian), e gives both the magnitude 
and direction of attitude error. In addition, for small perturbations. 


is a good approximation to the kinematic equation, and the state equation (54) 
becomes 
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e = to 


ea 


w = h 
ea 


There are many techniques for synthesizing the control law h. A simple 
example is 


h = - K,e - K 0 u (55) 

1 ^ ea 

where the constant gain matrices and are selected to provide the required 
bandwidth and damping in each axis. 

With the control law (55), the input to the moment trimmap becomes 

w T = A co - {if, e + |X, + 5(a) )]co } (56) 

al ao a 1 L 2 v oa J J ea v 

(shown schematically in fig. 20). The dynamic element (s) is included in 
the attitude perturbation controller to provide high-frequency cutoff. Thus, 
a feedback is closed around the process uncertainties, e, in equation (49) so 
that 


a 



(57) 


is sufficiently accurate if the angular acceleration command, co e , commanded by 
the attitude command generator is admissible, namely, if (A as ,co c ,co c ) is 
flyable and if the bandwidth of co e is suitably restricted. 


Now, consider the translation perturbation controller discussed at the 
end of the preceding section. The Coriolis effects may be included as follows. 
Let the matrix A define an axis system with respect to inertial space, and let 
E e = A{R SC - R s ) and^Fg = A(_V SQ - 7 S ) be the position and velocity errors, 
respectively. Then A = S(w)A, where w is the angular velocity of A, and 

R = S (m)i? + V 

e e e 


Hence, 


R e = 5(m)i? e + S(u>)[S(oj)l? e + V & ] + S(u)V e + 4(F gc - Fj 
K = Kc + aT ^ s ^ r 6 + S 2 t*)R e + 2S(oo)F e - R & ] 


The designer is free to choose A and R e . For example, let the error relax 
according to the linear law, 


S 

e 


-K,R - K„R 
l e 2 e 


Then the input to the trimmap is given by 
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Kl = Ka + + S 2 W + K 2 S(u) * K x ]R e + [25(c) + KJVJ 


In particular, if A = A vs , which aligns the first axis with the commanded 
velocity, V so> and maintains the second axis horizontal, then 


0) 


= -04 


vs 


* S l 4 3 )S( V f '.o /V 2 


where 


k = 


X T 

S J3%c 

■ (6 u ) : 

v 3 SQ J 


U 


SQ 




The results of simulation tests suggest that all significant coriolis effects 
are accounted for by the approximation in which k - 0 and the gains in 
figure 19 are replaced as follows, 

K 1 * + K 2 S ^V ) + K 1 

K 2 + 2 S(u v ) + K 2 


and K 3 is unchanged. 


Trajectory Command Generator 

The last two major blocks of the proposed structure of automatic flight- 
control systems are the trajectory command generator and the attitude command 
generator. Their function is to provide only admissible commands to the 
corresponding acceleration controllers. In this section, only the trajectory 
command generator is discussed. Since, within the scope of this report, the 
two generators may be considered to be very similar, the discussion applies 
also to the attitude command generator. 

In the hierarchy of control logic, the command generator is one level 
above the acceleration controller and one level below the generalized air 
traffic control (which includes the pilot) . Sufficiently smooth commands can 
be passed unmodified to the acceleration controller. However, in general, 
discontinuities will be present: the air traffic control may request a dis- 

continuous change in trajectory; the pilot may switch to a different control 
mode; the set of active sensors may change; or a strong disturbance due to 
wind or a partial failure may force the aircraft too far from the commanded 
trajectory to be brought back by the perturbation controller. In such cases, 
the command generator must generate an acceptable transition (flare) trajec- 
tory that returns the aircraft on target. The transition may be generated by 
means of a dynamical system as shown in figure 21. 
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Figure 21.- Generation of transition 
(flare) trajectory. 


The output of the command generator is 
given by 



where the quantities ( )* are the ATC 
command (see eq. (1)), and e p, ey, and 
ey are the modifications of that com- 
mand in position, velocity, and 
acceleration. 


Let the differential equation of the transition dynamics be 

e = h{e) 


with an output map 



(59) 


(60) 


If equation (59) is asymptotically stable and H{ 0) = 0, the output of the com- 
mand generator will approach the ATC command with time. To have continuity in 
commanded position, velocity, and acceleration, the dimension of e must be at 
least 9, that is, 3 for each axis. If initial conditions are chosen so that 


/e R {0)\ 


V°) 




/VO) 

I 

V s (°) 



(61) 


then, at the initiation of the transition, the command coincides with the 
actual position, velocity, and acceleration of the aircraft. 

The detailed shape of the transition is controlled by means of the func- 
tion h O) in the state equation (59). Generally, the state space will consist 
of at least two regions, one of which includes the origin e = 0. In this 
region, the function /z(e) may be linear. Thus, for example, let the small 
transitions be generated by three uncoupled, linear systems with constant 
coefficients , 
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( 62 ) 


where, for each i. = 1,2,3, the dimension of is 3 and the dimension of the 
constant matrices is 3 x 3, If the initial conditions are defined by the 
rows of the matrix 


e(0) = A VQ [e R {0) ,e v {0) (63) 

and the output map is defined by 

[e R (t) ,e v (t) = A^ Q e{t) (64) 

then the transition dynamics will be approximately invariant with respect to 
the commanded velocity axes given by the matrix A VQm Since the acceleration 
controller tracks the output of the command generator with small error, equa- 
tion (62) represents approximately the transition dynamics with respect to the 
longitudinal (i = 1), lateral (i = 2) , and normal (£ = 3) axes of the aircraft, 
respectively . The bandwidth of the transition can be made compatible with the 
restrictions of acceleration controller by a proper choice of matrices . 

Outside a neighborhood of e = 0, the function h{e) must be modified; 
otherwise, the magnitude restrictions of the acceleration controller will be 
violated. In this region of the state space of e , the design of 7z(e) may be 
based on such considerations as the optimization of transit time or transit 
energy with hard constraints on e. 

In effect, trajectory tracking errors have been sorted into three levels. 
Small errors are corrected by the perturbation controller without reinitializ- 
ing the command generator. Medium errors are corrected by means of the com- 
mand generator with linear transition dynamics. Large errors are corrected by 
means of the command generator with nonlinear dynamics. 

The total output of the command generator is given by equation (58). The 
generalized ATC command when not provided explicitly as a func- 

tion of time by ATC, must be generated onboard from a set of trajectory 
parameters that are either signaled by ATC or selected by the pilot. 

Finally, the command generator must contain a subblock within which 
autopilot modes can be defined for the control logic. The essential function 
of the mode variable is to specify which parameters of the commanded trajec- 
tory are to be tracked. A very simple example is given in table 2. 

TABLE 2.- EXAMPLE OF A MODE VARIABLE 


Trajectory parameter 

Axis 

to be tracked 

Longitudinal 

Lateral 

Vertical 

Acceleration 

0 

0 

0 

Velocity 

1 

1 

1 

Position 

2 

2 

2 


29 



The mode variable M in this case is three-dimensional. Each coordinate 
can take one of three values. Thus, there are 27 possible modes: M - (2,2,2) 
specifies position tracking in all three axes; M = (0,1,2) specifies the 
tracking of longitudinal acceleration, lateral velocity, and vertical posi- 
tion; and M = (1,1,1) represents velocity vector tracking mode, etc. Of 
course, other definitions of the mode variable are possible and are being 
investigated. Some of the commonly used modes can be included within the 
proposed structure by simply changing the set of active sensors. For example, 
if compass heading and barometric altitude are to be tracked instead of iner- 
tial heading and altitude, then compass and baro altimeter should be used for 
feedback instead of, say, MLS. Other modes, such as when the automatic flight- 
control system is allowed only limited authority and must interact with the 
pilot in the loop may be more difficult to include within the proposed struc- 
ture, but present indications are that such inclusions are possible. For the 
present purposes, however, it is sufficient to note that the automatic control 
logic must include a mode definition subblock. 

The proposed structure of the command generator is outlined in 
figure 22. 


Trajectory 

parameter 

select 


Mode 

select 


From sensors 



Figure 22.- Structure of the trajectory command generator. 
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PROPOSED STRUCTURE OF AUTOMATIC FLIGHT-CONTROL SYSTEM 


The overall logical structure of the automatic flight-control system 
developed here is outlined in figure 23. The structure consists of five major 
subsystems, namely, the trimmap, wind filter, attitude and throttle control 
systems, perturbation controller, and command generator. 



Figure 23.- Proposed structure of the automatic flight-control system 

for the AWJSRA. 

The decision to include a trimmap is motivated by the need to provide 
automatic envelope limiting and by the impracticality of overcoming the 
highly nonlinear characteristics of the aircraft by means of high-gain feed- 
back. In the trimmap, a priori information concerning the aircraft character- 
istics is used to generate open- loop control commands that trim the aircraft 
to a given acceleration vector. 

The decision to include a wind filter is dictated by the fact that aero- 
dynamic forces and moments are functions of the aircraft velocity relative to 
the air mass. 

Considerations of reliability and simplicity motivated the decision to 
impose a hierarchy on the control logic. The six degrees of freedom of the 
rigid aircraft are partitioned into a three-dimensional translation system and 
a three-dimensional rotational system. The function of the attitude control 
system is to execute the attitude commands provided by the (translation) 
trimmap. The bandwidth of the attitude control system is an order of magni- 
tude higher than the bandwidth of the translation control system. 
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Process uncertainties are controlled by means of a perturbation controller 
which closes the loop around the trimmap, wind filter, and attitude and 
throttle control systems. The design and implementation of the perturbation 
controller are drastically simplified by the decision to close the feedback 
through the trimmap. 

The subsystem composed of the perturbation controller, trimmap, attitude 
and throttle control system, and wind filter is an acceleration controller. 

Its input-output relation between the commanded acceleration V SQ and actual 
aircraft acceleration V s is approximately an identity everywhere on the flight 
envelope for suitably restricted acceleration commands. The function of the 
command generator is to give only admissible commands to the acceleration 
controller and to provide the interface between the control logic and the air 
traffic control or the pilot. 

As stated in the introduction, the purpose of the present report is not 
to present a complete design of an automatic flight control system, but, 
rather, to outline a structure of such systems. The discussion in the report 
leads to the structure composed of five major subsystems which are intercon- 
nected as indicated in figure 23. Some of the details within these subsystems 
discussed in the report are intended primarily to further clarify the purpose 
of each subsystem rather than as final designs. Indeed, the detailed struc- 
ture of each of the five subsystems is currently being developed and the 
results will be reported in forthcoming publications. However, the feasibility 
of the proposed structure has been tested by application to a simulation of the 
unmodified DHC-C8A and the Augmentor Wing Jet STOL Research Aircraft. The 
proposed logical structure has been shown to be feasible, and flight test 
evaluation will occur in the near future. 


CONCLUSIONS 


The proposed approach to the design of automatic flight control systems 
for advanced aircraft has several advantages, among which are the following. 

The approach is applicable to a large class of aircraft. 

The approach is nearly algorithmic. 

The tracking accuracy enters as an independent variable which may be 
varied over a wide range. 

There is an effective trade-off between tracking accuracy and flight com- 
puter requirements . 

Because the approach leads to a hierarchical system, questions of relia- 
bility are tractable. 
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The approach has been shown to be feasible, and flight test evaluation will 
occur in the near future. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, February 12, 1975 
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